On the physical basis for the nematic rubber elastic free energy 
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We discuss why it is physical to keep terms in the nematic rubber elastic free energy that reflect 
the order parameter dependence of the natural size of the network polymers. We address a point of 
difficulty in some mathematical approaches to this problem. 

PACS numbers: 61.30.Vx , 83.80.Va, 83.80.Xz, 62.20.Dc and 61.41.+e 



I. INTRODUCTION 



Nematic elastomers are formed by crosslinking nematic 
polymers. Such molecules order orientationally because 
of rod-like elements incorporated in their main chains 
(MC polymers) or pendant as side chains (SC). In both 
cases orientational order induces the backbones to change 
their mean shapes from spherical to elongated (prolate) 
or flattened (oblate) forms. Since chain shape and the 
macroscopic shape of a network are intimately related, 
there is a strong coupling between the nematic order 
and the shape of such solids. This obtains for both pro- 
late and oblate chain shapes - our conclusions will hold 
for both types of polymer. Several elastic phenomena, 
unique to nematic elastomers, arise from this coupling. 
The most relevant in a discussion of changing the mag- 
nitude of the order Q, is that such changes induce huge 
elongations and contractions (of many 100s%). The order 
can be changed by changing temperature, or by illumi- 
nation if at least some of the rods contain chromophores. 
The effect is reversible and is described in a monograph 
[H that summarises the work of many research groups. 
The magnitude of the spontaneous elongation, A,„, is pro- 
posed (as explained below) as an indicator of the shape 
anisotropy of the polymers, A™ — where f|| 

and £_L are the Flory effective step lengths parallel and 
perpendicular to the nematic ordering direction, n. 

Our purpose here is to clarify the way in which, at 
constant temperature, the nematic can perhaps change 
the magnitude as well as the direction of its order as dis- 
tortions are imposed. For instance, close to the nematic- 
isotropic transition, or in regions of extreme distortion 
close to the core of a disclination, Q may reduce. Some 
investigations, particularly from the mathematical com- 
munity, have had understandable concerns that the stan- 
dard nematic rubber elastic free energy apparently has 
terms unbounded from below in the limit Q ^ 1. To re- 
tain a physical description, authors have instead rigidly 
constrained such terms to be constant. We point out that 
(i) realistic chain models do have the critical terms vary- 
ing, and (ii) that there are more practical and physically- 
founded ways to deal with this apparent difficulty. 

A simple extension[l| to nematic elastomers of the clas- 
sical Gaussian theory of rubber elasticity yields a model. 
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where ji is the shear modulus of the rubber in the 
isotropic state. |o and In are the Flory effective step 

length tensors that give the mean square dimensions of 
a nematic Gaussian polymer, and there characterise the 
distribution of chain shapes. The former is at formation 
and the latter is that currently pertaining, that is after 
any director rotation or changes in Q induced by strain, 
or due to temperature or illumination change. Thus: 



{RM = le±L, (2) 



with £|| and £± the effective step lengths along and and 
perpendicular to the director, and where L is the chem- 
ical arc length of the polymer. £q and |„ depend on the 

order at formation and that currently pertaining (Q). 
The deformation gradient tensor is A and takes the body 

from its formation state to that current. Since rubber is 
a soft solid, deformations are strictly at constant volume 
and thus Det[A'^] = 1. Such a description of the net- 
work chains presumes they are long enough to be Gaus- 
sian, albeit anisotropic. It is reasonable to make this 
restriction for otherwise chains would not be dominated 
by their random configurations and nor would they be so 
highly extensible as their experimental response clearly 
demands. With this generic assumption, it then does not 
matter greatly what kind of model one adopts for chains 
- the Flory philosophy is that local molecular structure is 
washed out by randomness and can be encoded by £ (here 
two numbers £\\ and £±). The simplest model is that of 
the freely jointed chain with links of length a whereupon 
orientational order induces: 



Det[|] 



= (1 + 2Q), ^ = (1-Q), = (1 + 2Q)(1-Q)2 

(3) 

For elastomers over a large range of order parame- 
ters and hence also extensions the freely jointed chain 
model has proven highly accurate. It is found Q that 
spontaneous deformations arising from eqn ([l]). Am = 

1 /3 

((1 + 2Q)/{1 — Q)) correlates perfectly with indepen- 
dent measurements of Q from optical anisotropy which, 
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along with Am varies as temperature is changed. Here, 
distortions are with respect to a high temperature 
reference state where Q and A™ 1. Chains can 
be highly extended by high ncmatic order, attaining for 
instance a hairpin state which has a much more rapid 
increase of £|| with Q and hence with T In any event 
the picture of the effective step lengths and hence Det [|] 

varying with Q persists for all known chain models. It is 
unphysical to constrain the determinant to a fixed value. 
A supplementary argument that has been advanced for 
the fixing of Det[|] is that this is effectively proportional 

to the volume in spaced covered by the chain and that 
chains are incompressible, i.e. the Det[|] cannot change. 

The expression refers to the extent of the chain, a vol- 
ume proportional to N^^-^ since it is a random walk of 
lineal dimension proportional to n^/^. The actual vol- 
ume occupied by the monomers of the chain scales like 
N (times a^) and thus the density of segments of a par- 
ticular chain is ~ N~^/'^a~^ - chains are very dilute in 
a melt and most of the space they span is occupied by 
other chains. There is accordingly no volume constraint 
presented to a given chain. In fact this argument is at 
the heart of why chains adopt ideal statistics in a melt 
i- 

What then stops the free energy density ([1]) minimising 
at / = — oo by tending to Q = 1 where Dct[|] = and 

hence In ^Dct[|o]/Det[|]^ — ^ ~ool One must consider 

the dominant free energy in the problem, namely that 
of the underlying nematic phase rather than the weaker 
rubber part. The nematic chains have, in the absence of 
linking, a nematic free energy density fwiQ) that one can 
show to be additive to the rubber part arising on link- 
age. Model forms arise from simple Landau-de Genncs 
phenomenology, or from particular molecular models, for 
instance the extension of Maier-Saupc nematic theory to 
worm-like chains. In fact a freely jointed rod model of nc- 
matic chains would have exactly the Maier-Saupe ne- 
matic free energy since the rods are independent, except 
through their nematic interactions. As with the chains, 
it is not vital for our argument what model one adopts. 
They all have in common a free energy of order fcBTni per 
monomer where T^i is a characteristic temperature where 
nematic order is lost. The rubber free energy by contrast 
is of the order of k^T per network strand since the pref- 
actor in the free energy density eqn ([T]) is the modulus 
/i = UsksT where tis is the number density of strands. 
Thus the rubber component of the free energy density 
is of order 1/A^ smaller than the nematic part. N is 
the number of effective step lengths separating crosslink 

points along a chain. N has to be large, ~ 10 100 in 

order that chains are Gaussian and the network is rub- 



bery, that is highly extensible and dominated by entropy. 
The overall dependence of the free energy on Q is domi- 
nated by the minimum determining the stable state of the 
liquid under the same conditions. It is little perturbed 
by the Q-dependent terms arising from the rubber, in- 
cluding the lnDet[. . .Jterm. The rubber adds terms like 
Ji{~QqQ^ + Q^) [II, S for the ideal case and where 
spontaneous extension/contraction has been allowed to 
occur. (Qo is the order at crosslinking.) These terms 
cause minor shifts in the transition temperature and la- 
tent heat. In the non-ideal case, involving memory of 
orientation pertaining at formation, a term like —jfQ 
arises, which resembles the form of an external field and 
removes the isotropic high temperature phase Q identi- 
cally equal to zero altogether. 

Strictly speaking, even though the liquid ncmatic free 
energy dominates the nematic contributions from the 
rubber, the lnDet[. ..] term would still yield a sepa- 
rate minimum at Q = 1- However the use of the result 
£j_ = a{l — Q) near the point Q = 1 is also unphysical. 
This limit of perfect orientational order would suggest 
that (i?^) 0. Only for chains mathematically narrow 
and without any directional fluctuations (at T = 0) could 
attain this limit. In any case the Gaussian approximation 
of this result will also have long since failed. 

How should one proceed? Physically, one is con- 
strained to being near in Q to a modified nematic mini- 
mum in the free energy. One can thus reasonably ignore 
any minimum near Q = 1 and rely on the nematic free 
energy to naturally constrain Q-dependence arising in 
fc\- Practically, unless one is close to the N-I transition 
where the minimum in / near Q = is not stiff, one can 
take Q to be fixed in magnitude. If strains are extreme, 
then one can let Q have the freedom to respond to per- 
turbations arising from fci{X,Q) by taking the nematic 

free energy to be ^fq • {Q — Qmin)^ (with /q the curva- 
ture of the potential at Qmin), which was an approached 
followed in Q . These authors were concerned changes in 
the linear modulus from applied strains inducing changes 
in Q. This work also considered the role of A m inducing 

changes in the biaxial order if the principal stretch direc- 
tion was not aligned with Hq. Such changes too must be 
accounted for in the changing lnDet[. . . ] term that is the 
subject of this note. 

In summary we have argued physically as to why an- 
other source of Q-dependent free energy terms from rub- 
ber elasticity must not be constrained to take a constant 
value. In fact their form relates closely to well tested 
models of chain statistics and for the Gaussian theory of 
ncmatic rubber elasticity to hold, these terms must be 
retained. 
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